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We theoretically study the twofold and fourfold symmetric anisotropic magnetoresis¬ 
tance (AMR) effects of ferromagnets. We here use the two-current model for a system 
consisting of a conduction state and localized d states. The localized d states are ob¬ 
tained from a Hamiltonian with a spin-orbit interaction, an exchange field, and a crystal 
held. From the model, we first derive general expressions for the coefficient of the twofold 
symmetric term (Cf) and that of the fourfold symmetric term (C 4 ) in the AMR ratio. 
I 11 the case of a strong ferromagnet, the dominant term in C 2 is proportional to the 
difference in the partial densities of states (PDOSs) at the Fermi energy (E F ) between 
the de and dy states, and that in C 4 is proportional to the difference in the PDOSs 
at E f among the de states. Using the dominant terms, we next analyze the experi¬ 
mental results for Fe 4 N, in which \Cf\ and |C 4 | increase with decreasing temperature. 
The experimental results can be reproduced by assuming that the tetragonal distortion 
increases with decreasing temperature. 


1. Introduction 

The anisotropic magnetoresistance (AMR) effect is a phenomenon in which the elec¬ 
trical resistivity depends on the relative angle between the magnetization (M) direction 
and the electric current (I) direction (see Fig. I). 1-7 ) The AMR effect has been studied 
extensively both experimentally and theoretically since 1857, when it was discovered by 
W. Thomson.^ The AMR ratio, which is the efficiency of the effect, is generally defined 
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by 


Ap(0) p{<t>) ~ p± 


( 1 ) 


P P± 

with p_L=p(7r/2). Here, 0 is the relative angle between the thermal average of the spin 
(S) ( oc —M) and I, and p(0) is the resistivity at 0. 

Experimentally, Ap(0)/p has been measured for various ferromagnets such as Fe, 4 ) 
Co, 4 ) Ni, 4 ) Ni-based alloys, 2 ) and half-metallic ferromagnets. 8-15 ) In addition, the AMR 
ratios of many ferromagnets have been observed to be 

Ap(0) 


P 


= c 0 + c 2 cos 20, 


( 2 ) 


where c 2 is the coefficient of the twofold symmetric term and Co is chosen to be c 2 so as 
to satisfy Ap(7r/2)/p=0. 10,11,16-18 ) 

Theoretically, expressions for Ap(0)/p have often been derived by using electric 
transport theory based on the two-current model with s-d scattering. 2,3 > 5 > 6 > 19 > 20 ) The 
s-d scattering means that the conduction electron (denoted as s) is scattered by im¬ 
purities into the localized d states (denoted as d) with the exchange field and the 
spin-orbit interaction. As a representative study, Campbell, Fert, and Jaoul (CFJ) 2 ) 
derived an expression for Ap(0)/p for strong ferromagnets 21 ) such as Ni-based alloys 
[see Eq. (E-10)]. On the basis of the CFJ model, 2 ) Malozemoff obtained an expression 
for Ap(0)/p for weak ferromagnets 21 ) as well as strong ferromagnets. 5,6 ) In addition, 
we extended the CFJ model 2 ) and the Malozemoff model 5,6 ) to a more general model, 
which could systematically explain the experimental results of Ap(0)/p for various fer¬ 
romagnets including half-metallic ferromagnets. 19 ) 

We also derived the analytic expression for Ap(0)/p given by Eq. (2). 20 ) We then 
showed that the twofold symmetric feature of Eq. (2) could be intuitively explained 
by considering the d states distorted by the spin-orbit interaction. Note here that the 
crystal held of the d states was not taken into account in the derivation of Eq. (2). 
The d states were indexed by M— 0, ±1, and ±2, with M being the magnetic quantum 
number of the 3d states. Furthermore, the partial density of states (PDOS) of each d 
state at the Fermi energy (E F ) was assumed to be constant regardless of the d state. 

Recently, the AMR ratios of several ferromagnets 22-31 ) have been experimentally 
observed to be 


— ^— - = Co + C 2 cos 20 + C\ cos 40. (3) 

P 


Here, C 2 (CJ) is the coefficient of the twofold (fourfold) symmetric term, and Cq is 
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chosen to be C 2 — C 4 so as to satisfy Ap(7r/2)/p=0. For example, \C 2 \ and |C 4 | for Fe 4 N 
increase with decreasing temperature T as shown later in Fig. n. 22 ~ 26 ) The coefficients 
C 2 and C 4 were measured to be C 2 =—0.0343 and 6^=0.00556 at T=4 K. 22 ^ 

The set of Co, C' 2 , and C 4 , however, has seldom been derived within the frame¬ 
work of transport theory and has often been represented by phenomenological expres¬ 
sions. 2 '~ 29,32,33 ) We anticipate that expressions for C 0 , C 2 , and C 4 obtained by transport 
theory will play an important role in the analysis and understanding of the AMR ef¬ 
fect. We also predict that the fourfold symmetric term in Eq. (3) may appear under 
the crystal held of the d states, which was neglected in the previous models 19 - 1 [i.e., Eq. 
( 2 )]- 

In this paper, we obtained C 2 and C 4 by extending our model 19,2 °) to one with a 
crystal held. We hrst performed a numerical calculation of C 2 and C 4 for a strong ferro- 
magnet using the d states, which were obtained by applying the exact diagonalization 
method (EDM) to a Hamiltonian of the d states with a crystal held. The result revealed 
that C 4 appears under a crystal held of tetragonal symmetry, whereas it vanishes under 
a crystal held of cubic symmetry. We next derived general expressions for the resistivity, 
C' 2 , and C 4 for ferromagnets with the tetragonal held using the d states, which were 
obtained by applying hrst- and second-order perturbation theory (PT) to the Hamil¬ 
tonian. From the expressions, we obtained expressions for C 2 and C 4 for the strong 
ferromagnet with the tetragonal held. The result showed that C 2 cos 20 is related to the 
real part of the probability amplitudes of the specihc hybridized states and C 4 cos 40 is 
related to the probabilities of the specihc hybridized states. In addition, we performed 
a simple analysis of the experimental results of C 2 and C 4 for Fe 4 N using the dominant 
terms in C 2 and C 4 obtained by PT. The experimental results could be reproduced by 
assuming that the tetragonal distortion increases with decreasing T. 

The present paper is organized as follows: In Sec. 2, we obtain wave functions of 
the localized d states by applying hrst- and second-order PT to the Hamiltonian of 
the localized d states. Using the wave functions, we derive general expressions for the 
resistivity, C 2 , and C 4 for ferromagnets. In Sec. 3, we obtain expressions for C 2 and C 4 for 
a strong ferromagnet from the above-mentioned C 2 and C 4 . In addition, we perform the 
numerical calculation of C 2 and C 4 using the d states, which are obtained by applying 
the EDM to the Hamiltonian. We then compare C 2 and C 4 obtained by PT and the 
respective values obtained by the EDM. In Sec. 4, we analyze the experimental results of 
C 2 and C 4 for Fe 4 N. The conclusion is presented in Sec. 5. In Appendix A, we show the 
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matrix of the Hamiltonian. In Appendix B, we give the zero-order states of the d states, 
which are obtained by performing the unitary transformation on the perturbation term. 
In Appendix C, we describe the overlap integrals of the s—d scattering rate. In Appendix 
D, we give an expression for the s-d scattering rate. Section E shows that the present 
Ap(0)/p (= 26 * 2 ) coincides with our previous model 19 ' 20 ' 1 and the CFJ model 21 under 
appropriate conditions. 

2. Theory 

In this section, we obtain general expressions for the resistivity, C 2 , and C 4 in a 
model in which I flows in the x direction and (S) (oc— M) lies in the xy plane (see 
Fig. 1). We here use the two-current model with s-d scattering in which the conduction 
electron is scattered into the localized d states by nonmagnetic impurities. 2 > 3 > 5 > 6 > 19 , 20 ) 
The d states are obtained by applying PT to a Hamiltonian of the d states. We also 
explain the numerical calculation method for C 2 and C4, in which the d states are 
obtained by applying the EDM to the Hamiltonian. 


Z 



Fig. 1 . (Color online) Sketch of the sample geometry. The current I flows in the x direction, the 
thermal average of the spin (S) (oc—AT) lies in the xy plane, and </> is the relative angle between the 
I direction and the (S) direction. 


2.1 Hamiltonian 

We first present the Hamiltonian H of the localized d states of a single atom 19,341 
in a ferromagnet with a spin-orbit interaction, an exchange field, and a crystal field 
of tetragonal symmetry. This crystal field represents the case that distortion in the z 
direction is added to the crystal field of cubic symmetry. 35 ) Note that C 4 appears under 







a crystal field of tetragonal symmetry, whereas it vanishes under a crystal held of cubic 


symmetry, as will be described in Sec. 3.2. 

The Hamiltonian H is expressed as 


n = n 0 + v, 

(4) 

"Ho = "Hcubic — H ■ S, 

(5) 

h = Vso T htetra; 

( 6 ) 


with 

cubic £ E £ (\xy,xM))(xy,xM) I + \yz,xM))(v z ,xM)\ + \xz,xM))(xz,x<r 

<j=zL 

+E,(\x 2 - y 2 , X „(0 )><* 2 - y 2 , XM)\ + |3 - 2 - r 2 , xM))(^ 2 ~ r 2 , X.(0)|)] , 
H so — XL - S, 

htetra = £ M ®2,X<r(0))(z3,X<r(0)| + \yZ,xM))(yZ,xM)\) 

a=± 

+8 1 \2>z 2 — r 2 , Xct(4>)){^z 2 — r 2 , Xa(4>)\ , 

and 


S — (S x , Sy, S z ), 

( 10 ) 


( 11 ) 

H = H(coscf), sin</>, 0), 

( 12 ) 


where H > 0. Here, S is the spin angular momentum and L is the orbital angular 
momentum. The spin quantum number S and the azimuthal quantum number L are 
chosen to be S—1/2 and L— 2 . 19 ^ 

The above terms are explained as follows: The term "H C ubic represents the crystal held 
of cubic symmetry. The term —H S is the Zeeman interaction due to the exchange held 
of the ferromagnet H, where H(x—M and Hcx(S). The term 14 0 is the spin-orbit inter¬ 
action, where A is the spin-orbit coupling constant. The term Vt e tra is an additional term 
to reproduce the crystal held of tetragonal symmetry. The state | i, x CT (</>)) is expressed 
by lb Xo-(0)) = |f) |x<t( 0))- The state |i) is the orbital state, defined by | xy)—xyf(r), 
\yz)—yzf(r), \xz)=xzf{r), \x 2 -y 2 ) = \(x 2 -y 2 )f(r), and \3z 3 -r 2 )=^E(3z 2 - r 2 )f(r), 
with f(r) being the radial part of the 3d orbital, where r=\Jx 2 + y 2 + z 2 . The states 
| xy), | yz), and | xz) are referred to as de orbitals and \x 2 — y 2 ) and |3z 2 — r 2 ) are referred 
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to as dy orbitals. The quantity E e is the energy level of \xy) and E y is that of \x 2 — y 2 ). 
The quantity A is defined as A =E 1 — E e , S e is the energy difference between \xz) (or 
| yz)) and | xy), and <5 7 is that between |3 z 2 — r 2 ) and \x 2 — y 2 ) (see Fig. 2). The state 
|x<r(0)) (<x=+, —) is the spin state, i.e., 


\x + m = t> +11», 

(13) 

|x_«.)} = d=(-e-«1t> + U>), 

(14) 


which are eigenstates of — H ■ S. Here, |x+(0)) (|x_(0))) denotes the up spin state 
(down spin state) for the case that the quantization axis is chosen along the direction 
of (S). The state | t) (| I)) represents the up spin state (down spin state) for the case 
that the quantization axis is chosen along the z axis. 


Energy 

A 


l Sy 


3z 2 ~r 2 ,(jJ 

x 2 -y 2 ,a) 


0 


%S e 


xz,o), \yz,o) 
xy,o) 


Fig. 2. Energy levels of the 3d states in the crystal field of tetragonal symmetry. The energy levels 
are measured from E e . 


Regarding the parameters, we assume the relations A /H <C 1, |A|/A <C 1, S e /A <C 
1, and 5 7 /A <C 1, bearing a typical ferromagnet in mind. In particular, H, A, and |A| 
are roughly set to H~1 eV, A~0.1 eV, and |A|~0.01 e V. 36 ’ 37) 

2.2 Wave functions of localized d states 

To obtain the wave functions of the d states, we apply first- and second-order PT 
to % of Eq. (4). Here, "Ho of Eq. (5) is the unperturbed term, while V of Eq. (6) is 
the perturbed term. When the matrix of TL is represented in the basis set \xy, x CT (0)), 
| yz, x CT (0)), | xz, x<t{4>)), \x 2 - y 2 , X<r{4>)), an d |3^ 2 - r 2 , y CT (0)), the unperturbed system 
is degenerate (see Table A-l in Appendix A). We therefore use PT for the case that 
the unperturbed system is degenerate. 38,39 ^ First, the unitary transformation is per- 
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formed for the subspace with the basis set \xy, Xo-(<f>)), and \xz,Xa(4 > )) as 

mentioned in Appendix B. As a result, we obtain the zero-order states as |£ + , x + (</>)), 

14,X+(0)), |f-,X+(0)), |f+,X-(0)>> |4,X-(0)), and |f-,X-(0))- Here, f± and 5 £ repre¬ 
sent the eigenvalues of V in the above subspace, where is given by Eq. (B-l). The 
respective zero-order states are expressed as 

|£+,X+(0)) = A (6 e - y/5* + A 2 )| xy, x+(0)> + *A sin (f>\ yz, x+(</>)> - iX cos (j)\xz, x+(0)> 

(15) 

l^,X+(0)) = cos(%£,x+(0)) + sin (j)\xz, X+(0))> (16) 

|£-,X+(0)) = B (4 + V s * + X2 )\ x Vi X+(4>)) + *A sin <j>\yz, x+(<f>)) - i\ cos (f>\xz, X+{4>)) 

(17) 

lf+,X-(0)> = a (S £ - y/6% + A 2 )I xy,x-(</>)) -iAsin^l yz,x~(&)) + i A cos (j>\xz,x -(</>)) 

(18) 

14,X-(0)> = cos(/>|^,x_(0)) +sin0|x^,x_(0)), (19) 

|f-,X-(0)> = B (S e + v^T+ A 2 )|xr/, X-(0)) ^*Asin0|j/z,x-(0)) + iA cos 0|xz, X-(0)> 

( 20 ) 

with 


A = {251 + 2A 2 - 24+ A 2 ) -1/2 , (21) 

B = (2<5 2 + 2A 2 + 24V^s + A 2 )" 1/2 . (22) 

Next, using the basis set |f+,X±(0)), |4,X±(</>)), |f-,X±(0))> W ~ V 2 , X±(4>))> and 
|3^ 2 — r 2 , X±(0)), we construct the matrix of "H of Eq. (4) as shown in Table I. In the 
construction, we perform, for example, the following operations: 

^{L X S X + L y Sy)\^±,X-{4 > )) — — 7 

V / 24 2 + 2A 2 T 24 v ^ 2 TA^ 

X |^T _ I 3 ^ 2 _ r2 )( _icos2 ^lx+(^)) — sin 2(f>\x~{(f>))) 
+*y|® 2 -y 2 ,X+{<!>)) - ^A 2 | xy,x-{<t>)) 

+4&=f >/?+**) 

x[|x^)(f sin0|x+(0)) - cos0|x_(0))) 
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Table I. Matrix representation of H of Eq. (4) in the basis set |£ + ,x+(0)), 15 e , y+(</>)), |£-,X+(<^)), 
l£+>X-(0)}> I Se,X-(4>)), l£-,X-(0)}, \x 2 -y 2 ,X+(<P)), \x 2 -y 2 ,X-(0)), |3a 2 - r 2 ,x+(^)), and |3z 2 - 
r 2 , x-((/>))■ Here, £± is defined as £±=(<5 e ± yM 2 + A 2 )/2 [see Eq. (B-1)]. In addition, A and I? are given 


by Eqs. (21) and (22), respectively. In this table, {(f>) in y>(<(>) is omitted due to limited space. 



lf+.X+> 

l$e>X+) 

l£-.x+> 

l£+.x-> 

l4,X-> 

1C—> x— > 

2 2 
f - y , 

- id, 

|3z* - r z , 

|3z 2 - r 2 , 








x+> 

X— > 

X+> 

X-> 

<£+,x+l 

-§ + u 

0 

0 

0 

AA 

2 

0 

0 

—iXA 








x (2£_ — A) 



X(2£_ + A) 

X sin 20 

X cos 20 

<-5e,X+l 

0 

-f +«5e 

0 

AA 

2 

0 

XB 

2 

-zA 

1 2 

0 

a V3X 

1 2 

x/3A 

2 





x(-2^_ + 


x (—2£_ + 



X COS 20 

X sin 20 





A) 


A) 





(£->x+l 

0 

0 

-f+ £- 

0 

A B 

2 

0 

0 

-iXB 

V3X 2 B 

2 

-• 73X 2 B 

4 2 






X (2£+ — A) 



x (2£+ + A) 

X sin 20 

X cos 20 

(£+>x4 

0 

A A 

2 

0 

#+£+ 

0 

0 

—iXA 

0 

„• v'3>. 2 A 

4 2 

v/3A 2 A 

2 



X (—2^_ + 

A) 





X (2{_ + A) 


x cos 20 

X sin 20 

<5e,x-l 

A A 

2 

0 

AB 

2 

0 

f +*e 

0 

0 

;A 

2 

\/3A 

2 

■ \/3A 

L 2 


x(2£_-A) 


x (2£+ — A) 






X sin 20 

X cos 20 

<£_,X-I 

0 

A B 

2 

0 

0 

0 

f + £- 

-iXB 

0 

-• n/3X 2 B 
* 2 

V / 3A 2 B 

2 



X (— 2^_ + 

A) 





x (2£+ + A) 


x cos 20 

X sin 20 

{x 2 - y 2 , 

0 

1 2 

0 

iXA 

0 

iXB 

-f+A 

0 

0 

0 

x+l 




X (2£_ + |) 


X (2^+ + #) 





(x 2 - V 2 , 

iXA 

0 

iXB 

0 


0 

0 

f + A 

0 

0 

x-l 

X (2£_ + A) 


X (2£+ + A) 








(3z 2 - r 2 , 

\/3\ 2 A 

2 

a V3A 

1 2 

V3X 2 B 

2 

V3\ 2 A 

1 2 

V3X 

2 

-• Vzx 2 b 

L 2 

0 

0 

7 f+ A + 

0 

X+l 

X sin 20 

X COS 20 

X sin 20 

x cos 20 

X sin 20 

X cos 20 



(i'y 


(3 z 2 - r 2 , 

,• y/3 X 2 A 

1 2 

x/3A 

2 

■ V3X 2 b 

L 2 

x/3A 2 A 

2 

a \/3A 
* 2 

V3X 2 B 

2 

0 

0 

0 

■7^ + A + <5-y 

X-l 

x cos 20 

X sin 20 

X cos 20 

X sin 20 

X COS 20 

X sin 20 






Myz)(i cos0|x+(0)) + sin0|x-(0)))] j, (23) 

\(L X S X + LySy)\5 £ ,x~(</>)) = i^^\3z 2 - r 2 )(i sin 2(j)\x+((p)) - cos 20|x _{</>))) 

-i^\x 2 -y 2 ,x-((p)) + xy,x+((j)))- (24) 


Equations (23) and (24) play an important role in C 2 and C 4 as described in the (j) 
dependence of the wave functions in this section. 

Applying the usual Erst- and second-order PT to % in Table I, we obtain |i, X?(0))> 
where i (0 denotes the orbital index (spin index) of the dominant state in |i,x?(0))- 
The d state of the np spin | i,x+{<j>)) is expressed as 

lC+,X+(0)) = % +t+ (\x 2 -y 2 ,X-((/>))+wl+£ r2 + sm2(/)\3z 2 -r 2 ,x+((j>)) 


+wl+^ r 2 _ cos 2(j)\3z 2 - r 2 , x ~(</>))) • • • , 


(25) 
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%,X+{4>)) 


|£-,X+0)) 


x 2 -y 2 ,x+{4>)) 


|3 z 2 -r\ X+ m 


<A,+ ( !■ X 2 -y 2 ,X+ (</>)) + «&%,+ cos20|3^r 2 — r 2 , x+(^)) 
+wt?L r 2_sm2(j)\3z 2 -r 2 ,x-((j>)))--- , (26) 

c^_,+ (|^ 2 — y 2 , x_(0)) + w^~ 2 ^ r 2 ^ + sin 2(j)\3z 2 — r 2 , x+(0)) 
+ w 3z2_ r 2_ c °s2(/.| 3 z 2 -r 2 ,X-((/>))).- - , (27) 

Cx 2 -y*,+ (I X 2 - y 2 , X+(4>)) + W fz*-rC+ COS 2 ^l 3 ^ 2 “ ’ X+(0)» • • ■ t 

( 28 ) 

Qsz 2 —r 2 ,+ (13^ 2 - r 2 , x+(0)) + w 3 x fj y C+ cos201a; 2 - y 2 , x+(0)» 

(29) 


and the d state of the down spin | i, x-(0)) is expressed as 

l£+,X-(0)) = c ?+ _(|x 2 -t/ 2 ,x+( 0))+w|jl r2)+ cos20|3z 2 -r 2 ,x+(0)) 

+w|+r r2) _sin20|3^ 2 -r 2 ,X-(0)))... , (30) 

14, X-(</>)) = c s s ,-(\x 2 — y 2 1 X-{4>)) + w& ^-r 2 ,+ sin 2013z 2 — r 2 , x+(0)) 

+ w 3z 2 -r 2 -COs20|32; 2 -r 2 ,X-(0)))... , (31) 

lf-»X-(0)) = c ? _,_(|a; 2 -y 2 ,x+(0))+w|:r r2]+ cos20|3^ 2 -r 2 ,x+(0)) 

+ w 3 z^- r 2 ,- sin 20|3z 2 — r 2 , X-(0))) • • • , (32) 

k 2 ~y 2 ,X -(</>)) = dr 2 -y 2 ,-(k 2 -y 2 ,X-(0))+W3Xfccos20|32 2 -r 2 ,X-(0)))-- - , 

(33) 


|3- 2 -r 2 ,X-(0)) = c 3z 2_ r 2 t _(\3z 2 - r 2 ,x-((p)) + w 3 x 2_ v 7i cos20|x 2 - y 2 , X-(0))) 


® 2/ 


(34) 


In the right-hand side of Eqs. (25) — (34), we specify only the |3^ 2 — r 2 ,Xo-(0)) and 
|a; 2 — y 2 , X<y (0)) terms because these states contribute to the present transport in which 
I flows in the x direction (see Appendix C). The dominant states in |£ + ,x<r(0)), 
%,Xo{4>)), l£-,X<r(0)), \x 2 -y 2 ,Xc t(0)), and \3z 2 -r 2 ,Xa(0)) are respectively written 
as \t+,xM)), \^e,xM)), |f-,X<r(0)), k 2 -y 2 ,xM)), and |3^ 2 - r 2 , x<r(0)), although 
they are not shown in Eqs. (25) — (27) and (30) —(32). The other states, except for the 
dominant state in each |i,x<;(0))> represent the slightly hybridized states clue to the 
spin-orbit interaction. The quantity iu*A cos 20 or sin 20 represents the probability 
amplitude of |j, Xo-(0)) normalized by Cj j? . Here, iu*A is the coefficient of the cos 20 or 
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(a) 


(b) 


<j> 0 7t IA nil 37t / 4 n 5 tc/A / 2 1k! 4 2n 



tj> 0 7t IA 7t 1 2 3^/4 71 5^/4 3^/2 17TIA 271 



Fig. 3. (Color online) (a) Schematic illustration of the 0 dependences of \xz,X-{<t>)), cos0|3z 2 — 
r 2 , x+ (</>)), and sin 0|3^ 2 — r 2 , X-(0)} in ALyS'yjzz, X-(0)) °f Eh- (A-l). The upper part shows the top 
view (looking down along the z axis) of \xz,x~{<t>)) ■ The middle part shows cos0|3z 2 — r , x+(0)) in 
the xy plane. The lower part shows sin0|3z 2 — r 2 , %_(</>)} in the xy plane. Here, |xz,X-(0) is shown 
by the sky-blue or sky-blue-bordered orbital and |3z 2 — r 2 ,x±(0)} is shown by the yellow or yellow- 
bordered orbital. The blue curve in the middle part is cos 0 and that in the lower part is sin <f> . Note 
that the 0 dependent coefficients of |3z 2 — r 2 , X±(0)) are given by only cos 0 and sin0; that is, the 
prefactor of cos 0 or sin 0 is ignored. In addition, the color-filled orbitals (white orbitals with a colored 
border) express regions with a negative sign (positive sign) in the wave function, where the 0 dependent 
coefficients are taken into consideration in regard to |3z 2 — r 2 , X±(0))- (b) Schematic illustration of the 
0 dependences of |yz,X-(0)), sin0|3z 2 - r 2 ,x+(0)}, and cos0|3z 2 - r 2 ,X-(0)) in \L x S x \yz,X-{4>)) 
of Eq. (A-2). The upper part shows the top view (looking down along the z axis) of |yz,X-(</>)}■ The 
middle part shows sin0|3z 2 — r 2 , x+(0)} i n the xy plane. The lower part shows cos 0|3 z 2 — r 2 , X-(0)) 
in the xy plane. Here, \yz, X- (0)) is shown by the pink or pink-bordered orbital. The blue curve in the 
middle part is sin0 and that in the lower part is cos 0. The other notation is the same as in (a). 


sin 20 term normalized by c^, while is the coefficient of the constant term, which 
does not depend on 0. Such iu*A cos 20 and w l £ a sin 20 generate the twofold and fourfold 
symmetric terms of Ap(0)/p as described in Sec. 2.5. 
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(a) 


(b) 


(c) 


(j) 0 n / 4 k! 2 3^/4 n 5k IA 3k/2 1k!A 2k 





Fig. 4. (Color online) Schematic illustration of the 0 dependences of the dominant states in die states 
and the hybridized states in Eqs. (30) —(32). The dominant states in (a), (b), and (c) are |£+,x-(0)) 
of Eq. (18), |6 e ,X-(0)) of Eq. (19), and |£_,X-(0)) of Eq. (20), respectively. The hybridized states are 
represented by expressions with a probability amplitude of cos 20 or sin20, i.e., cos 20|3z 2 — r 2 , X±(0)) 
and sin 20|3z 2 — r 2 , x±(0))j where the prefactor of cos 20 or sin 20 is ignored. In each panel, the upper 
part shows the top view (looking down along the z axis) of the dominant state. In (a) and (c), the middle 
part shows cos20|3z 2 — r 2 ,x+(0)) hr the xy plane, and the lower part shows sin20|3z 2 — r 2 ,X-(0)) 
in the xy plane. In (b), the middle part shows sin20|3z 2 — r 2 ,x-(0)} in the xy plane, and the lower 
part shows cos 20|3z 2 — r, x+(0)) hr the xy plane. The state | yz, X-(0)} is shown by the pink or pink- 
bordered orbital, |a;z,X-(0)) is shown by the sky-blue or sky-blue-bordered orbital, and \xy, X-(0)) is 
shown by the gray or gray-bordered orbital. The state |3z 2 — r 2 , X-(0)) is represented by the yellow or 
yellow-bordered orbital. The color-filled orbitals (white orbitals with a colored border) express regions 
with a negative sign (positive sign) in the wave function including the probability amplitude. In the 
middle part of (a) and (c) and the lower part of (b), the blue and red curves are cos 20 and cos 2 20 
[=(1 + cos40)/2], respectively. In the lower parts of (a) and (c) and the middle part of (b), the blue 
and red curves are sin 20 and sin 2 20 [=(1 — cos 401/2], respectively. 
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Fig. 5. (Color online) Schematic illustration of the (f> dependences of the dominant states in dy 
states and the hybridized states in Eqs. (33) and (34). The dominant states in (a) and (b) are \x 2 — 
2 / 2 ,X— (0)) an d |3 z 2 — y 2 ,X-(f)), respectively. The hybridized states are represented by expressions 
with a probability amplitude of cos 2b, i.e., cos2b|3z 2 — r 2 ,X-(b)} and cos20|a; 2 — y 2 ,X-(b)}, where 
the prefactor of cos 20 is ignored. In each panel, the upper part shows the top view (looking down 
along the z axis) of the dominant state. In (a), the lower part shows cos20|3z 2 — r 2 ,X-(0)} in the xy 
plane. In (b), the lower part shows sm2c/)\x 2 — y 2 ,X-(0)) in the xy plane. The state |3z 2 — r 2 ,x'-(0)} 
is shown by the yellow or yellow-bordered orbital. The state \x 2 — y 2 , X-(</>)) is shown by the green or 
green-bordered orbital. The color-filled orbitals (white orbitals with a colored border) express regions 
with a negative sign (positive sign) in the wave function, where the probability amplitude is taken into 
consideration in regard to |3z 2 — r 2 ,X-(0)) and \x 2 — y 2 ,X-(0)}- In the lower parts of (a) and (b), the 
blue and red curves are cos 20 and cos 2 20 [=(1 + cos40)/2], respectively. 


2.3 Origin of cos 20 and sin 20 terms in d states 

We explain the origin of the cos20 and sin20 terms in Eqs. (25) —(34). In the 
de states, the cos 20 and sin 20 terms appear through de — d'y hybridization. In 
the d'y states, they appear owing to d'y — de — d'y 1 hybridization, in which the d'y 
states are hybridized to the d'y 1 states via the de states. These hybridizations are 
due to the specific matrix elements in Table I, i.e., (3 z 2 — r 2 ,yy(0)|'H|£±,Xcr(0)) 
and (3z 2 — r 2 , Xa'(4 > )\Td\S £ , Xo-(0)), with cr=+, — and cr'=+, —. We now focus 
on (3z 2 - r 2 ,Xa'{f)\'H\f±,X-((p)) and (3^ 2 - r 2 ,yv(0)|%|4,X-(0)>, where (3^ 2 - 
r2 iXa’iffdil^X+if)) and (3 z 2 — r 2 , X (r '(0)|'H|5 e , X+(0)) c& n also be discussed in a 
similar way. These matrix elements originate from only the cos 20 and sin 20 terms 
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in Eqs. (23) and (24). Here, cos 20 and sin 20 are formed by the multiplication of the 
following coefficients: 

(i) The 0 dependent coefficients of \3z 2 — r 2 ,X-(0)} in XL y S y \xz, x~(0)) of Eq. (A-l) 

and \L x S x \yz,X-(<j>)) of Eq. (A-2). 

Note that the 0 dependent coefficients of \x 2 — y 2 ,X-(0)} in Eqs. (A-l) and (A-2) 

are not responsible for the cos 20 and sin 20 terms in Eqs. (23) and (24). 

(ii) The 0 dependent coefficients of | xz,X-(4>)) and | yz,x~(<f>)) in Eqs. (18) —(20). 

We discuss (i). We first emphasize that the operations that generate \3z 2 — r 2 , Xo-(0)) are 
\L y Sy\xz, X-(0)) of Eq. (A-l) and \L x S x \yz, X-(0)) of Eq. (A-2). In Fig. 3(a), we show 
the 0 dependences of cos0|3^ 2 —r 2 , x+(0)) and sin 0|3z 2 —r 2 , X-(0)) in \L y S y \xz, X-(0)) 
of Eq. (A-l). Here, the coefficients of \3z 2 — r 2 , x±(0)) are given by only cos 0 and sin 0; 
that is, the prefactor of cos 0 or sin 0 is ignored for simplicity. When 0=0, the coefficient 
of \3z 2 — r 2 , X+(0)) is finite, whereas that of \3z 2 — r 2 , X-(0)) is zero. In brief, since the 
spin direction of | X— (0)) in \xz, X-(0)) is the x direction, S^lx-W) becomes S'Jx-W) 
=—||x+(0)). Namely, the spin is reversed by the operation of S y . In contrast, when 
0=7 t/ 2, the coefficient of |3z 2 —r 2 , x~(7t/ 2)) is finite, whereas that of \3z 2 —r 2 , x+(tt/ 2)) 
is zero. In short, since the spin direction of |x-(vr/2)) in | xz, X-(vr/2)) is the y direction, 
Sy\x-{Tt/2)) becomes S y |x-(tt/2)) =— ||x_(vr/2)). Namely, the spin is conserved under 
the operation of S y . In a similar way, we can consider the 0 dependence of the coefficient 
of |3^ 2 - r 2 , x±(0)) in >^L x S x \yz, X-(0)) of Eq. (A-2) [also see Fig. 3(b)], 

2.4 Illustration of d states 

In Figs. 4 and 5, we show schematic illustrations of the 0 dependences of the 
dominant states and hybridized states in Eqs. (30) —(34). The dominant states are 
l£+,X-(0)) of Eq. (18), 1 4,X-(0)) of Eq. (19), |f_,X-(0)> of Eq. (20), |x 2 - y 2 , X-(0)), 
and |3z 2 — r 2 , x~(0))- The hybridized states are represented by expressions with a prob¬ 
ability amplitude of cos 20 or sin 20, i.e., cos20|3^ 2 — r 2 , x±(0)), sin 20|3^ 2 — r 2 , x±(0)), 
and cos 201a; 2 — y 2 , X-(0)), where the prefactor of cos 20 or sin 20 is ignored for simplic¬ 
ity. Each probability is also given by cos 2 20 [=(l+cos40)/2] or sin 2 20 [=(1—cos 40)/2], 
Such 0 dependences originate from the 0 dependent coefficients of |3^ 2 — 'r 2 ,Xo-(0)) in 
\(L x S x + L y S y )\f±,X-(f)) of Eq. (23) and \(L X S X +L y S y )\S e , X-(0)) of Eq. (24). These 
operations are commented on as follows: 

(i) X(L x S x + L y S y )\f,±, X- (0)) Of Eq. (23) 

When 0=0, only A L y S y \xz, X-(0)) in this operation generates |3^ 2 — r 2 , X+(0)) [ see 
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the case of 0=0 in Figs. 4(a) and 4(c)]. This feature comes from the case of 0=0 
in Fig. 3(a). When 0= tt/2, only XL x S x \yz, x_(7t/2)) in this operation generates 
\3z 2 — r 2 ,x~ y(7 t/2)) [see the 0=7t/2 case in Figs. 4(a) and 4(c)]. This feature is due 
to the case of 0=7r/2 in Fig. 3(b). 

(ii) X(L x S x + L y Sy)\S e , X-(0)) of Eq. (24) 

When 0=0, only XL x S x \yz, X-(0)) in this operation generates |3 z 2 — r 2 ,X-(0)) [see 
the case of 0=0 in Fig. 4(b)]. This feature stems from the case of 0=0 in Fig. 
3(b). When 0=7t/2, only XL y S y \xz, x_(7t/2)) in this operation generates \3z 2 — 
r 2 ,y_( 7t/2)) [see the case of 0= 7t/2 in Fig. 4(b)]. This feature is due to the case of 
0=7r/2 in Fig. 3(a). 

Here, |au/,X-(0)) in |£±,X-(0)) and |<5 e ,X_(0)) is not responsible for the 0 dependent 
coefficients of |3z 2 —r 2 , x-(0)) as found from the fact that (3z 2 — r 2 , x o -'(0)|'H| xy, x<r(0)) 
does not depend on 0 (see Table A-l). 


2. 5 General expression for resistivity 

Using |i, x?(0)) of Eqs. (25) —(34), we can obtain a general expression for p(0). The 
resistivity p(0) is first described by the two-current model, 2 ' i.e., 

P+0)P-0) 


P(0) = 


(35) 


p+O) + p-0)' 

The quantity Po-(0) is the resistivity of the a spin at 0 with cr=+, —, where cr=+ (—) 
denotes the up spin (down spin) for the case in which the quantization axis is chosen 
along the direction of ( S ). The resistivity p c 00) is written as 


P<t(0) = 


m„ 


(36) 


n <T e 2 r <T (0)’ 

where e is the electric charge and n a (m*) is the number density (effective mass) of the 
electrons in the conduction band of the a spin. 40,41 ' The conduction band consists of 
the s, p, and conductive d states. 19 ' In addition, l/r CT (0) is the scattering rate of the 
conduction electron of the a spin, expressed as 


Ft(0) 


/s - cr ,• r=4- - 


?= + ,- 






with 


2vr 




(37) 


(38) 


Fs,(t—> dj,s(0) ^ 

where *=^+, S e , x 2 — y 2 , and 3^ 2 — r 2 . Here, 1 /r St(T is the s~s scattering rate, which 
is proportional to the PDOS of the conduction state of the o spin at E F , D < f\ 19 ' > The 
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s—s scattering means that the conduction electron of the o spin is scattered into the 
conduction state of the o spin by nonmagnetic impurities or phonons. The quantity 
1/T s ^^d i ,<;{ ( l ) ) i s the s-d scattering rate. 19,20 ) The s-d scattering represents the scattering 
of the conduction electron of the a spin into the a spin state in the localized d state of i 
and q by nonmagnetic impurities. The quantities i and q respectively denote the orbital 
and spin indexes of the dominant state in | i, x ? (0)). The localized d states | i, yy(0)) are 
given by Eqs. (25) —(34) obtained from RL of Eq. (4). The quantity represents the 
PDOS of the wave function of the tight-binding model for the d state of the i orbital 
and q spin at E F as was described in Ref. 19. 34 ) The conduction state of the a spin 
\e ikaX ,X(T (0)) is represented by the plane wave, i.e., \e lk,jX , x<r{4>))— (l/v / h2)e* fc<TX |y CT (0)), 
where k a is the Fermi wavevector of the a spin in the x direction (i.e., the I direction) 
and is the volume of the system. The quantitiy Vi mp (i? n ) is the scattering potential 
at R„ due to a single impurity, where R n is the distance between the impurity and the 
nearest-neighbor host atom. 19 ) The quantity N n is the number of nearest-neighbor host 
atoms around a single impurity, 19 ) n imp is the number density of impurities, and h is 
the Planck constant h divided by 2n. 

On the basis of |(i, x s (0)|e* fc<rX , Xa(4>))\ 2 described in Appendix C, we obtain 
i n Eq. (37) up to the second order of X/H, A/A, A /{H ± A), 5 t /H, 
5t/ A, or St/(H ± A), with t—e or 7. Details are given in Appendix D. 

Using these results, we obtain p a (4>) of Eq. (36) as 

Pa{4>) = P0,a + P‘1, a COS 20 + p 4 ,<x COS 40, (39) 

where po,<r is the constant term, which is independent of 0, p 2 , a is the coefficient of the 
cos 20 term, and p^ a is that of the cos 40 term. These quantities are specified by 

Po,a = Pol + Pol, (40) 

P2,a = P 2 I + P 2 I, (41) 

P4,cr = P±l, (42) 

where v of pu’l (u— 0, 2, 4 and u=0, 1, 2) denotes the order of X/H, A/A, X/(H ± A), 
5 t /H, 5 t /A, or 5 t /(H ± A), with t—e or 7. The quantities pul are obtained as 
(°) _ ,3 1 

P 0,± Ps,± T >i 2 — y 2 ,± T ^p s ,±^3z 2 —r 2 ,±, 

Po.± — ^ fix) 0 2 4 2 p s ,±-.$ + ,± + 3p S|± ^i e ,± + A 2 B 2 p s ,±_j^_,±) 


(43) 


J. Phys. Soc. Jpn. 


A 


4 \H- A 


A 2 id 2 _ + — ) Ps,±^t + ,T + gPs,±->-<5 e ,=F + -B" ( <^+ + — Ps,±-a£_, 
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(i) 3 A / A A 

P2,± = 


sy (a a t h ) p»,±^-^,±), 


(44) 

(45) 


(2) 3 

= 4 
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2 V A 
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16 A(A =f H) 

A / 3 | 1 A(<5 7 — S e ) 

<L | 4|2 A 2 


( XA 2 l + + XB 2 l _ — 2 A 2 d_l + — 2B L d + l_) >p St± ^ x 2_ y 2 ± 


+ 


(yA_) [-4 2 (y - «+)<*- + s 2 (y - «-m+] 


A 2 


16 A(A tH) 


(—AA 2 Z + — A B 2 l_ + 2 A 2 l + d_ + 2 B 2 l_d + ) >p s ,±-*-3z 2 -r2,±, (46) 


(2) = 3_(y 

p4 - ± 32 l A 


A 


32 \Ht A 


(—A’ A p s ,±->-£ + ,± + Ps,±-a<s £ ,± — A B p s ,±^_,±) 
j (A A Ps,±-a£+,=F — Ps,=t-x5 £ ,=F + ^ B p S) ±^£_ }Z f) 



Ps , d=—>-:r 2 — y 2 , zb 


As, zb—>• 3z 2 —r 2 , zb) 


(47) 


with 


m; 

P s - CT = 17^2 


n a e*T S)a 


rrr 


Ps&—>di,<; ~ 9 

-n 




(48) 

(49) 
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d±=S £ + A /2 ± V4 2 + A 2 , (50) 

l± = A - S E ± v^+A2, ( 51) 

where A, B, and are respectively given by Eqs. (21), (22), and (B-l) and A 2 (A 2 + 
B 2 ) =1 and A (A 2 d_ + B 2 d + )= 1/2 are used. Here, p s%a is the s-s resistivity and p so -_s.d ii? 
is the s-d resistivity. The s-d scattering rate l/r s>(T ^.d i)? is defined by 


—= '^n imp N n V imp (R a ) 2 \(3z 2 - r 2 ,xM)\e lk ° z ,xM ))\ 2 ^ 
27T 1 


(d) 

,T 


— Al n . TV -v 2 D^ 

, 6 imD J v n u rr 1 - y r 


(52) 


with 


^cr hj m p(/? n )p,j, (53) 

where is given by Eq. (C-5). The overlap integral (3 z 2 — r 2 , x<r(4>)\e lkaZ , Xa (0)) i n 
Eq. (52) can be calculated using Eq. (C-1). Note that Eq. (52) has been introduced 
to investigate the relation between the present result and the previous results 2,19 ) (see 
Appendix E). Equation (52) was used in the previous models. 2,19 ) 

On the basis of Eqs. (36) —(39) and (C-6) — (C-8) and Appendix D, the features of 
P 2 ,cr cos 20 and p^ a cos 40 in Eq. (39) are described as follows: 

(i) The resistivity p- 2 ^cos2(j) [see Eqs. (39) and (41)] is related to the real part of the 
probability amplitudes of |3z 2 — r 2 , x CT (0)) and \x 2 — y 2 , x CT (0)), which are given by 
Re[ 7773 / 2_ r 2 J cos 20 and Re[iu *2 2 J cos 20, respectively [see Eqs. (D-3) and (D-6)]. 
The quantity cos 20 in the probability amplitude is shown in Figs. 4 and 5. 

(ii) The resistivity p4 i(T cos40 [see Eqs. (39) and (42)] is related to the probabilities of 
\3z 2 — r 2 , Xa (0)) and \x 2 — y 2 , x<t( 0)), which are given by \w^ 2 _ r2 J 2 (l ± cos40)/2 
and K’ a ? _ I/ 2 <T | 2 (l + cos40)/2, respectively [see Eqs. (D-4) and (D-7)]. The quantity 
(1 ± cos40)/2 in the probability is shown in Figs. 4 and 5. 


2.6 General expressions for C 2 and C 4 

Using Eqs. (1), (35), and (39) —(42), we obtain a general expression for Ap(0)/p up 
to the second order of A/if, A/A, A/(if±A), 8 t /H , 8 t / A, or 8 t /(H± A), with t—e or 7 . 
The AMR ratio Ap(0 )/p is explicitly expressed by the form Ap(0)/p=U 0 + C 2 cos 20 + 
U 4 cos 40 , where Cq=C - 2 — C 4 . The coefficients C 2 and C 4 are written as 


65 = 


P-2,1 + P 2 !- 


( 0 ) 

Po,+ (1) 


„(°) 

doJl) 


(pSI + pZ ) 2 W- 


00 ) 


W p 2 .- + ~W P2 >+ + -( 0 ) 
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1 / n (0) n (0) 

1 / Po,+ ( 2 ) Po- ( 2 ) 

(0) I (0) d2,- + (0) P2,+ 


P0I + P00 V/C 


do/ 
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P 0 ,+P 4 ,- 


„(°) n (0) 
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P o°- (Po?+ + Po°-) Po°l (Poll + Po°-) 2 Po°Ip ( o°- 

1 pgj + P^ / Poll (1) Poll (1) \ (55) 

2 (pS+pD 2 Ur p " p& p "V’ (} 

where is given by Eqs. (43) —(47). Using Eqs. (54), (55), and (43) —(47), we can 
investigate C 2 and U 4 for various ferromagnets. Also note that Ap(0)/p (=2C 2 ) of the 
present model coincides with that of our previous model 19 ) and that of the CFJ model 2 ) 
under appropriate conditions (see Appendix E). 


2 .7 Calculation method of U 2 and C 4 by exact diagonalization method 

As a different approach from PT, we perform a numerical calculation of C 2 and 
U 4 using the d states, which are obtained by applying the EDM to H in Table I. 
The first purpose of this approach is to find the crystal held that leads to C 4 ^0. The 
second purpose is to check the validity of the results obtained by PT (see Sec. 3). The 
calculation in the EDM is as follows: 


(i) We numerically obtain | i, X 4 (0)) in Eq. (38) by applying the EDM to H in Table I. 

(ii) Utilizing the obtained |i,yy(0)) and Table C-l, we numerically calculate 

\(hxM)\ eik ° x ixM))\ 2 of E q- (C-6). 

(iii) Using the calculated (i,;u(0)|e VXe(0)) , we obtain a numerical value for 
A p(<p)jp of Eq. (1) with Eqs. (35)—(38). The numerical values of Ap(0)/p and 
Ap(tr/4)/p are represented by /o and /, T /r, respectively. 

(iv) When the AMR ratio is expressed as Eq. (3), we have 



(56) 

(57) 

(58) 


(59) 
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3. Application to Strong Ferromagnets 

On the basis of C 2 of Eq. (54) and C 4 of Eq. (55), we obtain expressions for C 2 
and C 4 for a strong ferromagnet with D^—O and The coefficients C 2 and U 4 

are compared with those obtained by the EDM. In addition, from the results of the 
EDM we find that C * 4 appears under a crystal field of tetragonal symmetry, whereas it 
vanishes under a crystal held of cubic symmetry. 


3.1 Expressions for C 2 and U 4 

Using Eqs. (43) —(47), (54), and (55), we obtain expressions for U 2 and C 4 for a sim¬ 
ple system with _=D^_ and D ^_ y2 _—D^ z \_ r2 _■ The relation D ^_ y2 _=D^ Z l _ r 2 _ 


gives 


= 0 . 


(60) 


where p 2 f is given by Eq. (45). In addition, in accordance with previous studies 42 ^ we 
assume n + —n _, m* + =m*_, and v + =v-, where v + —v- is satisfied by setting k + =k- in 
Eqs. (53) and (C-5). The expressions for C 2 and U 4 are then written as 
1 


Ch = 


( 0 ) ( 0 ) 

Po,+ ( 2 ) Po,- ( 2 ) 


Po,l + Pol \Pof- Po, 


X\ r 7 -r el 
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(61) 
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Here, we have 


A H + A 


(62) 
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(64) 

(65) 

( 66 ) 


^7 r x 2 -y 2 - ^3 z 2 —r 2 ,— i 


(67) 
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where 

r . _ _ Ps ^ di '~ , ( 68 ) 
Ps,+ 

with i=£+, S e , £_, x 2 — y 2 , and 3 z 2 — r 2 . The resistivity Ps^i- is given by Eqs. (49) and 
(52), where a in p s . a ^d, - is unspecified because the cr dependences of n a) m*, and k a are 
ignored as noted above. Furthermore, we note that r satisfies the relation 7y_ oc D^]_ 
[see Eq. (52)]. 

On the basis of (i) and (ii) of Sec. 2.5, the features of C2cos20 and C 4 cos40 are 
described as follows: 

(i) The term C 2 cos 20 is related to the real part of the probability amplitudes of 
\3z 2 — 'r 2 ,x CT (0)) and \x 2 — 2/ 2 ,Xcr(0))> which are given by Re[tc*’) 2 _ j , 2 J cos20 and 
Re[w* 2 _ y2 J cos20, respectively [see Eqs. (D-3) and (D-6)]. Concretely, C 2 cos20 
contains a single p^\ cos 20 in the numerator of each term in C 2 cos 20 [see Eq. 
(61)]. Here, p^ a cos2(j) is related to the real part of the probability amplitudes of 
\3z 2 — r 2 , x<r(0)) and \x 2 — y 2 , x CT (0)) as noted in (i) of Sec. 2.5. 

(ii) The term C 4 cos40 is related to the probabilities of \3z 2 — r 2 ,x<r(0)) and \x 2 — 
U 2 , X<y{4>)), which are given by |^2_ r , 2)(T | 2 (l±cos40)/2 and \w ^_ y2 J 2 (l + cos40)/2, 
respectively [see Eqs. (D-4) and (D-7)]. Concretely, C 4 cos40 contains a single 
pfl cos 40 in the numerator of each term in C 4 cos 40 [see Eq. (62)]. Here, pfl cos 40 
is related to the probabilities of |3z 2 — r 2 , Xo-(0)) and \x 2 — y 2 , Xcr(0)) as noted in (ii) 
of Sec. 2.5. Also, c 4 of Eq. (63) arises from high-order processes of d'y — de — oh/, in 
which the dy states are hybridized to the dy' states via the de states. Such processes 
reflect the fact that there are no off-diagonal matrix elements in the subspace of 
the dq states (see Table I). 

We next determine the effective value of the undefined parameter A/5 7 by comparing 
C 4 obtained by PT with that obtained by the EDM. We here put 


r £ 1 = a(i + y), 

(69) 

r e 2 = r £ , 

(70) 


where r) represents the difference between r e \/r e and r^/'W Figure 6 shows the |A|/5 7 
dependence of C 4 of Eqs. (62) and (59) for the systems with H —1 eV, A=0.1 eV, 
A=—0.01 eV, 37 ) r— 0, r 7 =0.01, 43 ^ r £ /r 7 =1, and ?/=0, 1, and 2. Here, r —0 and r 7 =0.01 
are set on the basis of those for Fe 4 N. 19 ^ The range of | A| /<5 7 is roughly assumed to 
be 0.5 < |A|/<5 7 < 1.5 by consideration of the above parameters and 5 7 /A <C 1. At 
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each rj, U 4 obtained by PT decreases with decreasing |A|/5 7 because of c 4 oc (A /S 7 ) 2 . In 
contrast, C 4 obtained by the EDM is nearly constant. In particular, when each de state 
has the same PDOS at E F (i.e., r)= 0), C 4 of Eq. (62) for PT becomes C 4 =c' 4 , whereas C 4 
of Eq. (59) for the EDM is evaluated to be C 4 ~0 independently of |A|/5 7 . In addition, 
the difference in C 4 between PT and the EDM decreases with decreasing |A|/5 7 . From 
these results, the effective value of |A|/<5 7 for PT is considered to be |A|/5 7 ~l/2. In 
other words, the present PT is unsuitable for application to systems with |A|/5 7 > 1. 

With regard to C 4 , from now on we focus on the dominant term with 
[\/(H± A)] 2 or (A/A) 2 under the condition |A|/<5 7 rsj 1/2. Namely, we neglect c ' 4 with 
(A/5 7 ) 2 [(A/A) — \/(H + A)] 2 , which corresponds to high-order processes. The domi¬ 
nant term in C 4 is thus expressed as 


C* = 


3 T £l - r e2 
32 1 + r + ry 


1 


r + r. 


7 


A 


H -A 


(r + r 7 ) 


(71) 


As seen from Fig. 6, C 4 of Eq. (71) agrees relatively well with that obtained by the 
EDM with |A|/A 7 =l/2. 

Furthermore, we extract the dominant terms from U 2 of Eq. (61) and C 4 of Eq. (71) 
taking into account the relation of typical ferromagnets, \A/H\ <C 1. The dominant 
terms are 


C 2 = 


3 / A 


r e 1 


A J (1 + r + r 7 )(r + r 7 ) ’ 


Ch = 


A 


r £ 1 - r £ 2 


32 \ Ay (1 + r + r 7 )(r + r 7 )' 




As a characteristic feature, U 2 of Eq. (72) is proportional to r 7 — r e4 (oc D\, _ — j s £l j, 


D 


(72) 

(73) 

(d) \ 


while C 4 of Eq. (73) is proportional to r £] — r e2 (oc D^f_ — D K ^_) 


)( d ) 


3.2 Various features of C 2 and C 4 

We investigate various features of C 2 and C 4 for a strong ferromagnet with H =1 eV 
and A=—0.01 eV. We here use C 2 of Eq. (61) and C 4 of Eq. (71) for PT and C 2 of Eq. 
(58) and C 4 of Eq. (59) for the EDM, where | A|/<5 £ =| A| /<5 7 =1 /2 is set for C 2 and C 4 
for the EDM. We also utilize Eqs. (69) and (70). As a particularly important result, we 
find that C 4 appears under the crystal held of tetragonal symmetry, whereas it vanishes 
under the crystal held of cubic symmetry. 44 ^ 

Using the EDM, we obtain the r £ /r 1 dependences of C 2 and C 4 for a system with 
the crystal held of cubic symmetry, where A=0.1 eV, 5 £ =^7=0, r— 0, r 7 =0.01, and ?/=0 
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Fig. 6. (Color online) The quantity |A|/<5 7 dependence of Ca for the systems with the crystal field 
of tetragonal symmetry. We here set H =1 eV, A=0.1 eV, A=—0.01 eV, r=0, r 7 =0.01, r e /r 7 =l, and 
77 = 0 , 1, and 2. The solid curves represent C 4 of Eq. (71) for PT. The dot-dashed curves represent C 4 
of Eq. (62) for PT. The dashed curves represent C 4 of Eq. (59) for the EDM. 




Fig. 7. (Color online) The quantity r E /r 7 dependences of C 2 and C 4 for the system with the crystal 
field of cubic symmetry. We here set H =1 eV, A=0.1 eV, A=—0.01 eV, <5 e =<5 7 =0, r=0, r 7 =0.01, and 
77 = 0 . The solid lines represent C 2 of Eq. (58) and C 4 of Eq. (59) for the EDM. 


(see Fig. 7). We find that C 2 can be expressed as a linear function of r £ /r 7 . The sign of 
C 2 changes in the vicinity of r £ /r 7 ~ 1 . Furthermore, C 4 takes a value of almost 0. 

Figure 8 shows the r £ /r 7 or i] dependences of C 2 and C 4 for a system with the 
crystal field of tetragonal symmetry, where A=0.1 eV, r=0, and 7]=0 and 1. From the 
results of PT, we find C 2 ~0 for the system with r £ i/r 7 =r £ (l + ? 7 )/r 7 =l and C 2 ^ 0 for 
that with r £ i/r 7 =r £ (l +i))/r 1 ^ 1. This feature mainly reflects Eq. (72). We also obtain 
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Fig. 8. (Color online) The quantity r e /r 7 dependences of C 2 and C 4 for the system with the crystal 
field of tetragonal symmetry. We here set H =1 eV, A=0.1 eV, A=—0.01 eV, r=0, r 7 =0.01, and 77=0 
and 1. The solid lines represent C 2 of Eq. (61) and C 4 of Eq. (71) for PT. The dashed lines represent 
C 2 of Eq. (58) and C 4 of Eq. (59) for the EDM, where |A|/<$ e =|A|/<5 7 =l/2. 


64=0 for the system with ? 7=0 and C 4 ^0 for that with rj^O because of C 4 oc r £ 1 — r £ 2 
(= r £ rj). The coefficients C 2 and C 4 obtained by PT qualitatively agree well with those 
obtained by the EDM. 

In Fig. 9, we show C 2 for systems with the crystal held of tetragonal symmetry, where 
r £ /r 7 =1, 1.2, and 2 and r]= 0. Here, C 4 for PT takes a value of 0 because of 77 = 0 , and 
| C 4 1 for the EDM is much smaller than | (P 2 1 • The upper panel shows the r 7 dependence 
of C 2 for systems with A=0.1 eV and r=0. The middle panel shows the r dependence of 
C 2 for systems with A=0.1 eV and r 7 =0.01. The lower panel shows the A dependence 
of C 2 for systems with r=0 and r 7 =0.01. In the upper panel, when r £ /r 1 = 1, C 2 for 
PT is close to that for the CFJ model, i.e., C2=(l/2)Ap(0)/p=(3/8)(A/i/) 2 (a — 1), 
with a=r £ =r 7 (see Appendix E.2). In the middle and lower panels, C 2 for PT takes 
a value of almost 0 in the case of r e /r 7 =l. The sign of C 2 for PT is negative in the 
case of r £ /r 1 = 1.2 or 2. In addition, IC 2 I for PT increases with decreasing r or A and 
with increasing r e /r 7 . These features mainly reflect Eq. (72). In all panels, C 2 for PT 
qualitatively agrees well with that for the EDM. 

Figure 10 shows C 2 and C 4 for systems with the crystal held of tetragonal symmetry, 
where r £ /r 1 —1 and p=1, 3, and 6. The upper panel shows the r 7 dependences of C 2 and 
C 4 for systems with A=0.1 eV and r=0. The middle panel shows the r dependences 
of C 2 and C 4 for systems with A=0.1 eV and r 7 =0.01. The lower panel shows the A 
dependences of C 2 and C 4 for systems with r=0 and r 7 =0.01. I 11 all panels, the sign of 







J. Phys. Soc. Jpn. 




Black curves: r^ry= 1 
Red curves: rjry=\.2 
Blue curves: 

Solid curves: PT 
Dashed curves: EDM 
• : C 2 of CFJ model 


10 


U 


- 0.001 


- 0.002 




Black curves: rjrf= 1 

j 

Red curves: r,Jrf= 1.2 

i 

Blue curves: 


Solid curves: PT 


Dashed curves: EDM 


-0.005 



0.15 


Black curves: rjr^= 1 
Red curves: rjr-f=\.2 
Blue curves: rjrf= 2 
Solid curves: PT 
Dashed curves: EDM 


0.25 0.3 


Fig. 9. (Color online) The coefficient C 2 for strong ferromagnets with the crystal field of tetragonal 
symmetry. We here set H =1 eV, A=—0.01 eV, r e /r 7 =1, 1.2, and 2, and 77 = 0 . Upper panel: The quantity 
r 7 dependence of C 2 for the systems with A=0.1 eV and r=0. Middle panel: The quantity r dependence 
of C 2 for the systems with A=0.1 eV and r 7 =0.01. Lower panel: The energy A dependence of C 2 for 
the systems with r=0 and r 7 =0.01. The solid curves represent C 2 of Eq. (61) for PT. The dashed 
curves represent C 2 of Eq. (58) for the EDM, where |A|/<5 e =|A|/5 7 =l/2. The dots represent C 2 of the 
CFJ model . 2 - 1 Note that C 4 is not shown because C 4 for PT is 0 and C 4 for the EDM is much smaller 
than |C* 21 - 


C 2 for PT is negative, while that of C 4 for PT is positive. In addition, \C 2 \ and C 4 for 
PT increase with decreasing r 7 , r, or A and with increasing rj. Such features are mainly 
due to Eqs. (72) and (73). The coefficients C 2 and C 4 for PT qualitatively agree well 
with those for the EDM. 
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Fig. 10. (Color online) The coefficients C 2 and C 4 for strong ferromagnets with the crystal field of 
tetragonal symmetry. We here set H= 1 eV, A=—0.01 eV, r e /r 7 =1, and rj=l, 3, and 6 . Upper panel: 
The quantity r 7 dependences of C 2 and C 4 for the systems with A=0.1 eV and r=0. Middle panel: 
The quantity r dependences of C 2 and C 4 for the systems with A=0.1 eV and r 7 =0.01. Lower panel: 
The energy A dependences of C 2 and C 4 for the systems with r= 0 and r 7 =0.01. The solid curves show 
C 2 of Eq. (61) and C 4 of Eq. (71) for PT. The dashed curves represent C 2 of Eq. (58) and C 4 of Eq. 
(59) for the EDM, where |A|/<5 7 =l/2 and \\\/5 e =\/2. 


4. Simple Analysis of C 2 and C 4 for Fe 4 N 

Utilizing the above results, we perform a simple analysis of the experimental re- 
sults 22) for the T dependences of C 2 and U 4 for an Fe 4 N 45 “ 47 ) film on a MgO(OOl) 
substrate, where I flows along Fe 4 N [100]. The experimental results clearly show the 
difference in the behaviors between the low-temperature range of 4 K < T < 35 K and 
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the high-temperature range of 35 K < T < 300 K (see circles in Fig. 11). Here, we re¬ 
gard Fe 4 N as a strong ferromagnet with H,| c ?~0. 19 ^ In addition, we mainly focus on the 
effect of the PDOSs of the de states on C 4 . Note that we do not take into account the 
realistic crystal structure of Fe 4 N (i.e., a perovskite-type structure 45 -*) for simplicity. 48 ) 

From Eqs. (72) and (73), we first obtain simple expressions for C 2 and C 4 for Fe 4 N. 
By taking into account the relation for Fe 4 N, i.e., r < 1 and r 7 <C l, 19 ) C 2 and C 4 are 
given by 


where 


c 2 = - ± 


c 4 = 


K = 


3 ( A 
^A 

3 (\ 
32 VA 

3 ( A 
A 


r e 1 k(1 - Re i) k(1 - Re 2 - A R e ) 


r + r~ 


rei - r e2 


1 + R 

- 


p, 


■e2j 


R= — = 


' r + r 1 

1 + R 

2 

77) Dl 

) n el ~ -> 

77) A2 

-n-e 2 — - 

r 7 

r 7 

Ps,—/Ps,+ 

Ps- 


1 + R 

1 + R ’ 


Ps^d-y,— / Ps,-\- Ps^d-y,— 


(74) 

(75) 


(76) 


with Here, is proportional to 

which is the difference in the PDOSs at E F among the de states. 

We next determine parameter sets for A, A, R , R ei , R e2 , and A/L that can reproduce 
the experimental result for the T dependences of C 2 and C 4 . The quantity A is set to 
A=—0.013 eV for Fe. 36 ) The quantity A is assumed to be A=0.1 eV. 37 ) Here, the T 
dependence of A is considered to be negligibly small, because the decrease in the lattice 
constant due to a decrease in T is less than 0.5%, 22 ) where A of 0.1 eV is due to the 
Coulomb interaction between a magnetic ion and the surrounding ions. We accordingly 
adopt the T dependences of R , R £4 , R e2 , and A R e . The T dependence of R (R e4 , R e2 , 
and A R e ) is shown in the middle (lower) panel of Fig. 11. Details of the parameter sets 
are given below. 

We express R of Eq. (76) as 



P s — ydry , — 


p(" p +rft 

imp 

P S —^ d'y , — 


(77) 


where p s - and p s ^ t d 1 - are assumed to be p s -=p 1 J™ p + pf/_ and p s _> c q,-=p!™f ( q - ( see 
Sec. 2.5). The quantity p),™ 13 is the s-s resistivity due to impurities and p!™?)/, _ is the 
s-d resistivity due to impurities, where d 7 , — represents the dy states of the down spin. 
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Here, p™ p / P™-^d n ,~ i s se t so that 


fillip 

n im P 

rs—>d-y,— 


« 1, 


(78) 


by considering that p'^ d ^ _ (p™ P ) satisfies the relation p^ d ^ _ oc D^_ (p]™ 0 oc D ( l' I49 ' > ) 
and also Fe 4 N satisfies D^l_ D^\ On the other hand, p ph _ is the s-s resistivity due to 
the phonons. This p^ h _ depends on T through the influence of the number of phonons, 
which depends on T. 

The parameter sets in the high- and low-temperature ranges are noted below. 


(i) In the high-temperature range of T > 35 K, we have 

R = OT, (79) 

A R s = 0, (80) 

C (T) 

R e2 = R fA = 1 - (1 + R*), (81) 

K 

with 0=0.0270, T *=35 K, R*=QT*, and '=—0.0126, where C A ' is the exper¬ 
imental value of C 2 at T=T*. 22 ) 


The procedure for determining this parameter set is as follows: First, C 4 is ex¬ 
perimentally observed to be almost 0. Since C 4 oc A R e [see Eq. (75)], we assume 
AR £ =0 (or R £ i=R e2 )] that is, the PDOSs of the de states at E F take the same value. 
From the viewpoint of the crystal structure of Fe 4 N, this assumption may imply 
that the crystal exhibits cubic symmetry. 48 ) Next, since |C* 2 | gradually decreases 
with increasing T, we straightforwardly take into consideration the T dependence 
of R of Eq. (77), where R is included in the denominator of C 2 of Eq. (74). The 
denominator 1 + R is expressed as 


1 + R pa 


1 + 


ph 

_fR_ 

imp ' 
' s — yd'y ,— 


(82) 


by using Eqs. (77) and (78). Here, pf'_ is assumed to be proportional to T on the 
basis of the experimental result for the T dependence of the total resistivity. 50 ) 
Thereby, R (« pt-fpsZp,-) is S iven b Y 


R = 9T, 


(83) 


where 9 is a constant number. On the other hand, R e2 is determined so that Eq. 
(74) satisfies the condition (T, C 2 )—(T*, •'). Namely, R £2 is expressed as _R £ 2=1 — 



J. Phys. Soc. Jpn. 


(C^ V K )(1 + OT*). Substituting this R e2 and AR £ =0 into Eq. (74), we obtain 


Co = 


cf" ) (l + 07’*) 

1 + 9T ' 


(84) 


From the least-square fitting of Eq. (84) to the experimental result for C 2 , we 
determine 6 to be 0 (=0.0270), 51 ^ where the fitting is done for T=35—150 K by 
paying attention to the relatively low temperature side. Using 0=0.0270, we can 
also evaluate R e2 of Eq. (81). 


(ii) In the low-temperatnre range of T < 35 K, we have 


R = 0T, 

(85) 

AD _ T -T* /1 + R ( Ti) \ 

ARe 7] - T* V k/4 ) ’ 

(86) 

p _ K e2 ~ K s 2 rp , K e2 1 ~ K e2 1 1 

P'* _ p 1 p'* _ p 1 ’ 

(87) 

Re 1 = Re2 + A R £l 

(88) 


with 7)=4 K, R* £ . 2 =1-{1 + R*)C ( P/k, R^l-il + R^C^+AC^ 0 )/^ Ri=6T h 
C^=— 0.0343, and C^ l) =0.00556, where ( C | T ^) is the experimental value of 

C 2 at T=Ti (C 4 at T=71). 22 ) 

The procedure for determining this parameter set is as follows: We first adopt 
R=QT , which is the same as Eq. (79) in the high-temperature range, on the basis 
of the expeimental result of the T dependence of the total resistivity. 50 Second, 
since C 4 was experimentally observed to be a linear function of T, we assume C 4 to 
be C 4 —pT + q, where p and q are constants. The constants p and q are determined 
so that Eq. (75) satisfies the condition (T, C 4 ) = (Ti, C^), (T*,0). As a result, C 4 
is expressed as C 4 ={T — T*)C^ r, ' > /(7) — T*). From this C 4 and Eq. (75), we obtain 
A R e of Eq. (86). The obtained AR e may indicate the following two properties: One 
is that the crystal has tetragonal symmetry, which generates ATLy^O due to the 
difference of the PDOS at E F among the de states. The other is that the tetragonal 
distortion increases with decreasing T. Third, R e2 is assumed to be R e2 =p'T + q' as 
a simple form, where p' and q' are constants. The constants p' and q' are determined 
so that Eq. (74) satisfies the condition (T, C 2 )—(T h C^), (T*,C 2 T J ). 

Substituting the above-mentioned A and A, Eqs. (79) —(81), and Eqs. (85) —(87) into 
Eqs. (74) and (75), we obtain C 2 and C 4 for Fe 4 N, where C 4 in the low-temperature 
range was described above. In the upper panel of Fig. 11, we show the T dependences 
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of C -2 and C 4 . We find that C 2 and C 4 for PT successfully reproduce the experimental 
results. In particular, the experimental results in the range 4 K < T < 35 K, in which 
the change of |C 2 | is about four times as large as that of |C 4 |, can be explained by the 
ratio of the coefficients of A R e between Eqs. (74) and (75). 

Finally, we comment on the above-mentioned T dependence of A R £ , i.e., the dif¬ 
ference in the PDOSs at E F among the de states. The T dependence of A R e has been 
assumed to arise from the increase of the tetragonal distortion due to a decrease in 
T. The tetragonal distortion may originate from the anisotropic thermal compression 
of the lattice. This compression is considered to be due to the adhesion between the 
Fe 4 N him and the MgO substrate. We expect that such an assumption will be verified 
experimentally in the future. 

5. Conclusions 

We theoretically studied the twofold and fourfold symmetric AMR effects of fer- 
romagnets. In particular, we obtained the coefficients of the twofold symmetric term 
(cos 20 term) and the fourfold symmetric term (cos 40 term) in the AMR ratio, denoted 
as C 2 and C 4 , respectively. We used the two-current model for the system consisting of 
the conduction state and localized d states. The localized d states were obtained from 
the Hamiltonian with the spin-orbit interaction, the exchange held, and the crystal 
held. Details are given as follows: 

(i) We performed the numerical calculation of C 2 and C 4 for a strong ferromagnet using 
d states, which were obtained by applying the EDM to the Hamiltonian. The result 
revealed that C 4 appears under the crystal held of tetragonal symmetry, whereas it 
vanishes under the crystal held of cubic symmetry. 

(ii) We derived general expressions for the resistivity, C 2 , and C 4 for ferromagnets with 
the tetragonal held using the d states, which were obtained by applying hrst- and 
second-order PT to the Hamiltonian. From the expressions, we obtained expressions 
for C 2 and C 4 for the strong ferromagnet with the tetragonal held. The result showed 
that C 2 cos 20 is related to the real part of the probability amplitudes of the specihc 
hybridized states \3z 2 — r 2 , x CT (0)) and \x 2 — y 2 , Xo-(0)) and C 4 cos 40 is related to the 
probabilities of |3 z 2 — r 2 ,Xcr(0)) and \x 2 — y 2 ,Xo-(0))- In addition, we investigated 
various features of C 2 and C 4 obtained by PT and found that they qualitatively 
agreed well with those obtained by the EDM. 

(iii) We analyzed the experimental results of the T dependences of C 2 and C 4 for an 
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Fe 4 N film on a MgO substrate using the dominant terms in C '2 and C 4 obtained by 
PT. The dominant term in C 2 was proportional to the difference in the PDOSs at 
E f between the de and dy states, and that in C 4 was proportional to the difference 
in the PDOSs at E F among the de states. The experimental results in the high- 
temperature range (35 K < T < 300 K) were well reproduced by taking into account 
the T dependence of the s-s resistivity and by assuming that the PDOSs of the de 
states at E F took the same value. This assumption might imply that the crystal 
structure of Fe 4 N exhibits cubic symmetry. Also, the experimental results in the low- 
temperature range (4 K < T < 35 K) were successfully reproduced by assuming that 
the difference in the PDOSs at E F among the de states increased with decreasing T. 
This assumption suggested that the tetragonal distortion increases with decreasing 
T. Here, the tetragonal distortion was considered to originate from the anisotropic 
thermal compression of the lattice due to the adhesion between the MgO substrate 
and Fe 4 N film. 
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Appendix A: Matrix Representation of EL 

We construct the matrix of EL of Eq. (4) as shown in Table A-l. 

In the construction, we perform, for example, the following operations: 


\L y S y \xz,X-{4>)) 


\Ly\xz)Sy\X-(4 >)) 

i\(-\x 2 -y 2 ) + V3\3z 2 - r 2 ))^(icos (j)\x+((j))) + sin </>\x -(</>))), 

(A-l) 


\L x \yz)S x \x-{4>)) 

i\{\x 2 - y 2 ) + V3\3z 2 - r 2 ))^(isin0|x+(0)) - cos<£| X-{<t>)))- 

(A-2) 


\L x S x \yz,X-{4>)) 
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Table A-l. Matrix representation of H of Eq. (4) in the basis set \xy, x+ (</>))> \yz, x+(<!>))> 

I xz,x+{4>)), \xy,x-{<t>D, \yz,x-{<t>)), \xz,x~(<t>)), \x 2 - y 2 , X+0)}, |3z 2 - r 2 , x+ (<!>)), \x 2 - y 2 , X-(<!>)), 
and |3 z 2 — r 2 ,x~ (</>))■ In this table, {(ft) in x<z{(ft) is omitted due to limited space. 



\xy,x+) 

1 yz,x+) 

\*z,x+) 

\xy >x-> 

1 yz, x- ) 

\ xz , X—> 

1®* - y\ 
x+) 

- r*, 

X+> 

P 2 - y 2 , 

X—> 

\ 3 z z - r z , 
X-) 

(xy,x+ 1 

H 

2 

i ^ sin 0 

—i -A cos 0 

0 

^ COS 0 

^ sin 0 

0 

0 

—iX 

0 

(yz,x+\ 

—i-j sin 0 

-f 

0 

— ^ COS 0 

0 

_lA 

1 2 

— 2 COS 0 

— 2 COS 0 

^ sin 0 

X'l A sin 0 

<*2,X+I 

i ^ cos 0 

0 

-%+Se 

— ^ sin 0 

;A 

1 2 

0 

sin 0 

2 ^ X sin 0 

— COS 0 

' / J* X COS 0 

(xy,x- 1 

0 

— COS 0 

— ^ sin 0 

H 

2 

-2-^ sin 0 

i ^ cos 0 

—iX 

0 

0 

0 

{yz, x-1 

COS 0 

0 

~ih. 

1 2 

i ^ sin 0 

f +^e 

0 

— ^ sin 0 

_ V|A sin< ^ 

2-5- COS 0 

2 COS 0 

{xz, X- 1 

^ sin 0 

1 2 

0 

cos 0 

0 

f + < 5 e 

COS 0 

— ^ x COS 0 

2-j sin 0 

—2 X^ X sin 0 

(x - y , 
x+l 

0 

i ^ cos 0 

2-j sin 0 

iX 

— sin 0 

COS 0 

-f+ A 

0 

0 

0 

( 3 z 2 - r 2 , 
X+l 

0 

i cos 0 

—i ^ x sin 0 

0 

— y/ ^ x sin 0 

— V ^ X COS 0 

0 

-f+A + 

0 

0 

(x 2 - y 2 , 

X- 

iX 

-j sin 0 

— -J COS 0 

0 

cos 0 

—i-j sin 0 

0 

0 

f+ A 

0 

— 

1 

0 

sin 0 

^ X COS 0 

0 

—i v ^ x cos 0 

i v ^ x sin 0 

0 

0 

0 

■^ + A + ( 5 -y 


Equations (A-l) and (A-2) play an important role in C 2 and C4, as described when we 
discuss the (ft dependence of the wave functions (see Sec. 2.2). 


Appendix B: Zero-Order States 

Performing the unitary transformation for V of Eq. (6), we obtain the zero-order 
states. 38 ^ 

Table B-l shows the matrix representation of V in the subspace with the basis 
set \xy, x<r{4 > ))i | yz,Xa(<ft)}, and \xz,Xa((ft)), with a=+ or —. The eigenvalues of V are 
obtained as £+, S £f and £_, with 


e± 


8 s ±y/5fn? 

2 


(BT) 


In the case of a—+, the eigenstates for £+, S e , and are respectively given by 
|£+,X+(0)) of Eq- (15), | S £ ,x+((ft)) Of Eq. (16), and |f-,X-(0)> of Eq. (17). In the 
case of a =—, the eigenstates for £ + , S £i and £_ are |£ + ,X-(0)) of Eq. (18), | S £ ,x~(<ft)) 
of Eq. (19), and |£_,x~(0)) of Eq. (20), respectively. These states correspond to the 
zero-order states in PT. 


Appendix C: Overlap Integral of s—d Scattering Rate 

We briefly discuss |(i, X<;(4 > )\ etkaX ^ X<r(0))| 2 in Eq. (38), where | i, x ? (</>)) is represented 
by a linear combination of | xy,Xa{<ft)), \v z , X<r{<ft)), \xz,Xa{4>)), W - y 2 ,Xa{<ft)), and 
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Table B-l. Matrix representation of V of Eq. (6) in the subspace with the basis set \xy, x± (</>)), 

I yz,x±(<j>)), and \xz,x±{4>)). _ 



\xy,x±((f>)) 

|y^,x±(0)) 

kux±(</0) 

(xy, x±(</>) 1 

0 

sin 0 

cos 0 

{yz,x±{4>) 1 

=F i\ sin 0 

4 

0 

{xz,x±{(t>)\ 

cos 0 

0 

<5, 


|3 z 2 - r 2 ,xM))- 

On the basis of a previous study, 3 ^ we first give the following overlap integral: 


(t*v, X<r'(0)|e ife<T ' r , xM)) = 


POO /»oo POO 


' —oo J —oo J —oo 


32vrrC 


/ (r) iiu—j=e l ^" T ' r dxdydzSo 

Ouu --IV On 


(C-l) 


ka ' r iXM)) = (l/V^)e lka ' r Xa((p) (see Sec. 2.5), k a =(k Xt<T ,k yta ,k ZtlT ), 


Vn(k* + c 2 ) 3 V ^ ki + c 

i _ 

Here, we have \e 

k(j ~~||, and |//i/, Xa'(<!>))=f(r)fJ’VX<T'i with n=x, y, z, u=x, y, z, cr=+, -, and cr'=+, 
—, where /(r) is the radial part of the 3d orbital expressed by f (r)=Te~ < ’ r , and T and 
C are constants. The state | fi>v,X<r , (<f>)) denotes the fus orbital with a' spin. 

Using Eq. (C-l), we can calculate the overlap integrals for realistic orbitals. In the 
case of k a —(k cr , 0,0), corresponding to I//x (see Sec. 2.5), we have 


(x 2 ~y\Xa' (</>)!' e lkaX ,Xa (<t>)) = 


(3z - r ,Xa'((/>)\e l ° x iXM)) = 


2y/?> 


da^o 


(C-2) 

(C-3) 


{xy,Xa'{(l))\e l aX ,xM)} = (yZiXa'i&le 1 ” x ,xM)) = (xz,Xa>{<t>) |e* ^iXM)) = 0, 


(C-4) 


with 


1927tTC^ 

Vn(K + cY 


(C-5) 


Equations (C-2) —(C-4) mean that only \2>z 2 — r 2 , x<j{4>)) and | x 2 — y 2 , Xa(4>)) contribute 
to the transport of I//x. 

Next, using Eqs. (C-2) —(C-5), we obtain | (i, Xv(<^)|e ifccr:E , Xa(<f>)) f i n Eq. (38). Here, 
|*,X?(<£)) is S iven simply by \i, x ? (0))=£E a TM)\ji xM), where a}£(0) is the 
coefficient of |J, X<r(0))- In the case of j=x 2 — y 2 or 3 z 2 — r 2 , aU(0) corresponds 
to q i? , Ci^Wj’1 cos 2(f>, or sin 20, as seen from Eqs. (25) —(34). As a result, 
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\{hxM)\ eikaX xM))\ is expressed as 


2 a x 2 -y 2 ,<r( ( l ) )9(r ^ ®32 2 —r 2 ,cr 


{<t>)9o 


(C-6) 


i • 12 

In addition, |(i, X^(0)|e ifccra; , Xo-(0)) leads to the following two types of expressions (see 
Table C-l). Type 1 is written as 

|4 + 4 cos20| = |4| 2 + gKcl 2 + (44 + 44) cos 20 + g I4h cos40, (C-7) 
where £ 0 (4) is the coefficient of the constant term (cos20 term). Type 2 is 


|4sin201 2 = i|4| 2 - i|4| 2 cos40, 


(C-8) 


where 4 is the coefficient of the sin 20 term. Type 1 generates the twofold and fourfold 
symmetric terms and type 2 generates the fourfold symmetric term. 


Appendix D: s—d Scattering Rate 

Using Eqs. (38), (25) —(34), and (C-2) — (C-8), we obtain the sum of the s-d scattering 
rates of the second term in the right-hand side of Eq. (37), i.e., JT ^2 l/r SjCr ^. ( i ii? (0). 
We can express JT ^ l/r S)CT ^ dii? (0) as 
1 


EE 




Ts,<r->-di,?(0) ^ 


2 7T 2 

= ~T _ ^'impA n Vimp(7?n) (^-0,ff T A 2<j>,a T A40 jCr ). 


(D-l) 


Here, X 0r7 is the constant term, which is independent of 0, X 2 ^ a is proportional to 
cos 20, and A 4 ^ j(T is proportional to cos 40. The terms of a=+ are as follows: 


1 12 2 

^+5 + 

| C E,+ | °2 

W 3 2 2_ r 2 + 


D 


(d) 

hr 


W 


1 I 2 2 

£— ,+ 

| C E,+ | °2 

W 3z 2_ r 2 t+ 


+ l C <5 e ,+ |^ ( °1 + 2°2 


+ \Cx 2 -y 2 ,+ V ( o\ + ~o\ 


+ I c 32 2 -r 2 ,+ | ( °2 + 2°1 


32 2 -r 2 ,+ 
2 


D 


(<0 

^e,+ 


D 


(<0 

€-,+ 


w. 


3z 2 —r 2 ,+ 
it 3z 2 —r 2 ,+ 

4 2 — y 2 ,+ 


(d) 


x 2 —y 2 ,+ 
(d) 

3z 2 —r 2 ,+ 


+ | C E,-| (°l + 9°: 


W. 


1 

2 


12 2 

S e ,~ 

1 °2 

W 3 2 2 -r 2 ,+ 


E - 

3z 2 —r 2 ,+ 
2 


D 


D 


(d) 

X- 


2 / 2 I ■ 
+ C£_- o ± + — o: 


w 


E>- 

3z 2 —r 2 ,+ 


D 


D 

D 

{d) 

Z+- 

(d) 

E,-> 


(D-2) 
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Table CM. Expressions for | (i, xR</>)|e lfe ^ x ,x<r(0))| 2 - Types 1 and 2 are given by Eqs. (C-7) and 
(C-8), respectively. Here, Eqs. (C-2)—(C-5) are used. 


S e 

Z- 


x 2 — y 2 


9 Or 


{hX<(4>)\e aX ,xM)} 


+ + 9% 

+ + g 2 + 
+ + g\ 
+ + g 2 + 


3z 2 -r 2 + + g\ 


£+ 

Z- 

z+ 

£- 

a: 2 — y 2 

3z 2 - r 2 

Z+ 

4 

Z- 


+ 

+ 

+ 


9 2 - 


9 2 - 


9- 


9 2 - 


9 2 - 


9 2 - 


9 2 - 


9 2 - 


+ 9i 
+ 9l 
+ 9+ 


5^ C £ + ,+™£ 2 -r 2 ,+ Sin2 <t> 


°Se,+ l I - 273^fc- r 2,+ COS 

2j3 C t-,+ W b-r*,+ 8in2< f> 


I i 1 x 2 —y 2 

C x 2_ y 2 + ( - - — W. 


1 


v - 2y/3 W 3z*-r*,+ COS 

f I 1 Q^ 2 _t* 2 -U 

C 32 2_ r2 , + + 2 W x 2_y2X COS 

(|-i73 w fe-^,- cos2< ^) 

^5 c s s ,+ w t 2+ -r 2 ,- s in20| 

C «-+ (H i^fe-r 2 ,- C°S 2 </>) 

^ c ?+,- w &:r2,-Sin20 


C <5e,- ( 9 - 

1 


1 


Se- 


2 VS^-^+COS 
2 


2?3 C *-,-™fa-r 2 ,- sin 20 


Cx 2 -y 2 - ( 2 


/1 1 x 2 —y 2 ,— 

\2 2\/3^32 2 —r 2 ,— COS 


C 32 2_ r 2 _ (-^5 + 


1 3 z 2 —r 2 — 


2 2 ’ COS 


2 w x 2 -y 2 , 


C t+<- ( I - 


2 2 ? 3 «&’-r 2 , + cos 20) 

. 12 
1 , ...Se,~ 


i73 c <5e,-w 3 £ ; 2 _ r 2 )+ Sin 20 

(l - 


2 i7gw- 2 ’_ r , 2|+ cos20) 


type 


2 

1 

2 

1 

1 

1 

2 

1 

2 

1 

2 

1 

1 

1 

2 

1 


^2<f>,+ — 20iC>2 


|c« 5 e ,+ | 2 Re 
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Df\ 

0 £ ,+ 


T" Ipr 2 —?/ 2 ,+ l Re 
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3 z 2 —r 2 ,+ 

| 32 2 —r 2 ,+ 
^ x 2 —y 2 ,+ 


D 


D 


( d ) 

X 2—y2 _\_ 

(d) 


32 2 —r 2 ,+ 


+ |c ?+ _| 2 Re 


w 


£+,- 

32 2 —r 2 ,+ 




+ |c^_ _ 2 Re 


w 


€- - 
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cos 20, 
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W. 
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|2 2 

„s e ,+ 

1 o 2 

W 3z 2 —r 2 ,+ 


D 


(d) 


I 2 2 


°2 

W 3z 2 -r 2 ,+ 


D 


(d) 
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|2 2 

a: 2 —j/ 2 ,+ 

°2 
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D 
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2 .2 
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D 


(d) 


s e ,~ 


w 


- 

3 z 2 -r 2 ,+ 


D 


(d) 


cos 40, 


(D- 4 ) 


with 0!=(a; 2 —y 2 , X(j(0)|e* fcCTa; , x<j(0)) and o 2 =(3z 2 —r 2 , x CT (0)|e* fc<Ta: , x CT (0)), where o x and 
o 2 are calculated in Eqs. (C-2) and (C-3), respectively. The terms of cr=— are as follows: 


*0,- — c £+,+ ( °i + 


w 


+ 2 l c<5 -^+l a 


2 2 
2 


w 


5e,+ 

3z 2 —r 2 ,— 


£+,+ 

3z 2 —r 2 ,— 

2 


D 


(d) 

(+,- 


D 


(d) 

5e,+ 


+ | C «-,+ | (°l + 2°2 


w 


+ o K+ - 


1 2 2 




t+ - 

3 z 2 -r 2 - 


*-,+ 

3z 2 -r 2 ,- 

2 


D 


(d) 

«-,+ 




+ l c 5e,-T ( °1 + 2°: 


w. 


+ 2 


2 2 
°2 


W 


3z 2 —r 2 ,— 


8e- 

3z 2 —r 2 ,— 
2 


D 


(d) 

5e,— 


D 


(d) 

- 


+ \Cx 2 —y 2 ,— | 2 ( df + -O 2 


w. 


3 z 2 -r 2 - 


D 


(d) 


+ I c 3 z 2 -r 2 ,-d ( °2 + 2°1 


3 z 2 -r 2 - 

ID 2 2 

x 


D 


(d) 

3z 2 —r 2 , — ’ 


(D- 5 ) 


X2 <f>~ — 2 01O2 


c $+,+ 2 Re 


w. 


£+,+ 

3 z 2 —r 2 ,— 


D 


(d) 

£+,+ 


+ |c^_ j+ 1 Re 
+ |c< 5 e _| 2 Re 


w. 


t-,+ 

3 z 2 -r 2 - 


D (d) 
^£-,+ 


w 


8e~ 

3 z 2 -r 2 - 


D 


(d) 

- 


T |c a: 2_.y2 j _| Re 


w. 


x 2 -y 2 ,- 
3 z 2 —r 2 - 


D 


(d) 



J. Phys. Soc. Jpn. 


+ |C3^- r 2 ? -| 2 Re ItfS-S'- UfA,- COs2 4>, 


(D-6) 


Xaa, - — 


m. I A 2 L,£+>+ V r>( d ) 


_ \ r , | 2 n 2 „,«+>+ nW 

“ 9 l C W+l °2 w 3z 2_ r 2_ ^ +i+ 


|„ |2 2 nM 

l C A+l °2 W 3z 2_ r 2_ ^Se,+ 
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~KA 2 °l\ W tC r >,-\ Dfl_ 
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+ |c ^-^,-| 2 02 I ^-y 2 ,- 
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cos 40. 


(D-7) 


Appendix E: Relation between the Present Model and Previous Models 

E. 1 Correspondence to our previous model 

We show that Ap(0)/p (= 26 * 2 ) of the present model coincides with that of our 
previous model 19 ) under the condition D^fl=D^\ which indicates that the orbital i 
dependence of the PDOS is ignored. Under this condition, p s . a -Hi, : q is replaced by p Si 0 —>d,<r 
[see Eqs. (49) and (52)]. This replacement leads to p^±=0 [see Eq. (45)]. 

On the basis of Eq. (7) in Ref. 19, we first give an expression for the resistivity with 
spin-flip scattering p CT ( 0 ), i.e., 

/ \ = P+WuW + P+0)P-,+ 0) + P-Q)P+,-Q) rF 1 \ 

PW) P + (0)+P-(0) + (l + «)P + .-(0) + (l + «- 1 )P-. + (0)’ 1 

with a—m*_n + /{m* + n-)^ where p + _( 0 ) [p_ j+ ( 0 )] is the resistivity of the spin-flip 
scattering from the up spin to the down spin (from the down spin to the up spin). 

Using Eqs. (1), (E-l), and (39) —(47), we can obtain an expression for Ap(0)/p. The 
coefficient C 2 is finally obtained as 


^2 = T (pg|W+P^2), 


A — [(Ps,+ + p S) H— >d,+ )(Ps- + Ps, -»d,- + P-,+) + {pa,- + Pa, ->d,-)P+,-. 
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x[Ps,+ + Ps,-\ — >d,+ + ps- + ps, — >d- + (1 + a)p + - + (1 + a 1 )p_ i+ ], (E-3) 

Yl = (ps- + Ps--td,-)(p8,- + Ps, - >d,~ + P-,+ - P+,-) 

+ (Ps,- + Ps, - >d,- + P-,+) [(1 + a )P+- + (1 + a 1 )P-,+]) (E-4) 

Yl = (p s ,+ + Ps,+^d,+ )(.Ps,+ + p s ,+-+d,+ + P+,- — P-,+) 

+ (Ps,+ + Ps,+^d,+ + P+,— ) [(1 + a )p+,- + (1 + a 1 )P—,+] 5 


~(2) 3 

P2+ = o 


~(2) 3 

P2- = o 


A 


H - A 


(Ps,-\ — >d,~ ps+ —»d,+); 


( Ps ,- >d,+ Ps, - +d,— )• 


(E-5) 
(E- 6 ) 

^ i vrs ’— rs ’—(E-7) 

Here, we have ignored the 0 dependences of p + ._(0) and p_. + (0) in the same manner 
as in Ref. 19; that is, p + _(0)=p+_ and p_ i+ (0)+p_ i+ have been used. 

On the assumption of H 3> A, we express O 2 of Eq. (E-2) as 

3 / A A 2 1 

0*2 = ~ J — >d,— Ps,-\ —>d,+)Al + ( Ps, >d,+ Ps, »d,—)U 2 ], (E-8) 

where [A /(H ± A)] 2 ~(A /H) 2 has been used. As a result, Ap(0)/p (=2C 2 ) becomes 
Ap(0) 3 / A A 2 1 


p 4 \ iif / P s H— >d,+)Yi 4“ {ps, - >d,+ Ps ,- >d— )^]- 


Equation (E-9) corresponds to Eq. (28) in Ref. 19. 


(E-9) 


E.2 Correspondence to CFJ model 

We show that Ap(0)/p (— 2 C 2 ) of the present model coincides with that of the CFJ 
model 2 ) under the condition of the CFJ model, i.e., D^— 0, r £ i=r £ 2 —r 7 =cq r <C a, and 
r <C I. 19 -* Here, C 2 is given by Eq. (61). 

Under the above condition and the assumption of H 3 > A, we express C 2 of Eq. 
(61) as 


c2 -jKu) (a_i) ’ 

where [\/(H ± A)] 2 ~(A/J/) 2 . As a result, Ap(0)/p (= 26 * 2 ) becomes 

Ap(0) 3 / A 


Equation (E-ll) is Ap(0)/p of the CFJ model in Ref. 2. 


(E-10) 


(E-ll) 
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Fig. 11. (Color online) Upper panel: The T dependences of C 2 and C 4 for Fe 4 N. The solid curves 
represent Eqs. (74) and (75) for PT. The dots represent the experimental values in the temperature 
range from 4 to 300 K for the case of J//Fe 4 N [100 ]. 22 ' > Middle panel: The T dependence of R in Eq. 
(76). The expression for R is given by Eqs. (79) and (85). Lower panel: The T dependences of R el , 
R e 2 , and A R e in Eq. (76). The expressions for R e 1 , R e2 , and A R e are given by Eqs. (81) and ( 88 ), 
(81) and (87), and (80) and ( 86 ), respectively. The black solid curve (black dashed curve) represents 
Rci (Re 2 ). The blue solid curve represents A R e . 
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48) We report a basic feature of Fe 4 N. The crystal structure of Fe 4 N is a perovskite- 
type structure, in which N is located at the body center position of fcc-Fe. 45 ^ The 
unit cell with a cubic shape consists of a corner site and three face-center sites. 
Here, the de states for each spin at the corner site are considered to be degenerate. 
The three face-center sites are in the xy, yz, and xz planes and are denoted F xy , 
F yz , and F xz , respectively. We can then specify the ground and excited states of 
the de orbitals at F xy , F yz , and F xz taking into account the effect of N at the body 
center site. At F xy , the ground state is the xy orbital, while the excited states are 
the yz and xz orbitals. At F yz , the ground state is the yz orbital, while the excited 
states are the xy and xz orbitals. At F xz , the ground state is the xz orbital, while 
the excited states are the xy and yz orbitals. In this system, the PDOSs of the xy, 
yz, and xz states for each spin at E F take the same value. In contrast, when the 
system has tetragonal distortion in the z direction, the PDOS of the xy state at 
E f is different from that of the yz or xz state at E F . 

49) We obtain p 1 s mp oc from Eqs. (17) and (19) in Ref. 19. 

50) From Fig. 2 in Ref. 22, we roughly evaluate the total resistivity p to tai to be 
Ptotai-P/3 + 14 /iff • cm. This p to tai is expressed as p to tai= p+p-/(p+ + P-) ~ p - 
owing to the relation for Fe 4 N, i.e., p~/p+ <C l. 19 ' Here, p_ is simply given by 
p~ = ps - + Ps—td,—— p s - + Ps- + Ps-td -■ bmce p s _ depends on 1 , we assume 
p^ h _—T/3 and p™^ + p n ^ d _ = 14. Namely, p v ^_ is proportional to T. 

51) In this study, we only consider the relation R oc T [see Eq. (83)] on the basis of 
Ref. 50. We here do not judge the validity of 0=0.0270. The present model, which 
consists of the d states of a single atom, does not take into account the d states 
in the unit cell of the realistic crystal structure (i.e., perovskite-type stucture). In 
such a model, it is inconsequential to judge the validity of the numerical value of 
0 . 



